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Abstract 



We employ an effective Lagrangian approach to derive the leading-logarithm two-loop 
electroweak contributions to the muon anomalous magnetic moment, a^. We show that 
these corrections can be obtained using known results on the anomalous dimensions of 
composite operators. We confirm the result of Czarnecki et al. for the bosonic part and 
present the complete sin 2 6\y dependence of the fermionic contribution. The approach is 
then used to compute the leading-logarithm three-loop electroweak contribution to 
Finally we derive, in a fairly model-independent way, the QED improvement of new- 
physics contributions to and to the electric dipole moment (EDM) of the electron. We 
find that the QED corrections reduce the effect of new physics at the electroweak scale 
by 6% (for a M ) and by 11% (for the electron EDM). 
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1 Introduction 



Historically, the measurement of the anomalous magnetic moment of the muon |I[ 

G „ = ^i^ = — ■ ^ — - - 1 = (11 659 230 ±84) x 1(T 10 (1) 

has provided the most convincing test of the validity of QED [0. As the precision in the 
experimental measurement and in the theoretical determination are progressively improving, 
the muon anomalous magnetic moment a M may soon offer a test of the Standard Model (SM) 
electroweak theory and of some of its possible extensions. 

From the experimental point of view, the E821 experiment at Brookhaven National Labo- 
ratory is expected to improve the accuracy in the a M measurement to the level of 4 x 10~ 10 , and 
possibly to 1-2 xl0~ 10 if large statistics is accumulated ||. Let us consider now the present 
status of the theoretical prediction for a M . It is convenient to separate the total result into 
several different parts. The pure QED contribution, which is known to order a 5 [[J, is 

fl QED = 5g4 7Q6 ± 2 ) x 10 -11 _ ( 2 ) 

The part affected by strong-interaction contributions, which contains the largest source of uncer- 
tainty, comes from the hadron vacuum polarization and the hadronic light-by-light amplitude. 
A recent analysis relating, by means of dispertion relations, the hadronic vacuum polarization 
to data from e + e~ annihilation and r decays gives 0] 

aj ad (vac pol) = (6951 ± 75) x 1CT 11 . (3) 

The error is dominated by the experimental uncertainty and will be significantly improved 
by future measurements at BEPC in Beijing, at DA$NE in Frascati, and at VEPP-2M in 
Novosibirsk. The effect of higher-order hadronic contributions have also been evaluated || 

aJ ad (h.o. vac pol) = (-101 ± 6) x 10~ u . (4) 

The most recent theoretical estimate of the hadronic light-by-light contribution gives 

aj ad ( 7 x 7 ) = (-79 ±15) x 10~ n . (5) 

Finally, the SM electroweak contribution at one loop is [|7| 

1 



a™(lloop) = ^^ 
V lJ 24V2V 2 



1 + -(1-As 2 w ) 2 



195 x lO" 11 , (6) 



where is the Fermi constant and s 2 ^ = sin 2 9w = 1 — Myy/Mg. As first noticed in 
ref. |§, two-loop electroweak corrections are quite substantial, because of large contributions 
0(G M m^a/nha (M/rrif)). Here M represents the W or Z mass and rrif indicates a light fermion 
mass. Combining together the contribution containing closed fermionic loops |], [K| and the 
one from the other relevant two- loop diagrams (usually indicated as bosonic part), Czarnecki 



et al. obtain 11 



a^ w (2 loop) = (-44 ± 4) x 1CT 11 , (7) 

for an Higgs mass Mg = 250 GeV, where the error is associated to the uncertainties in Mh, 
the hadronic contributions, and higher-loop effects. 

In this paper, we show how the leading 0(G^ m 2 a/ir In (M/rrif)) two- loop electroweak 
corrections to the anomalous magnetic moment, a^ w (2 1oop)ll, can be easily obtained with the 
help of effective theories and Wilson renormalization group. The relevant anomalous dimensions 
of composite operators can be extracted from known results of the QCD corrections to flavour- 
violating bottom-quark transitions mediated by the magnetic dipole operator JHJ. In this 
way, we derive a^ w (2 loop)i,L without directly computing any Feynman diagram and confirm 
the results of refs. ||, |10|, Moreover, since the renormalization-group technique actually 
includes all leading QED logarithms, we are able to give an analytic expression for the leading 
three-loop contribution to a^ w . This contribution turns out to be very small, but its knowledge 
allows to eliminate the uncertainty on aj^ w from higher-order effects. Including the different 
terms, we obtain 

a™ = (153 ±3) x 1CT 11 , (8) 



where the central value corresponds to Mh = 150 GeV. This leads to a theoretical prediction 
for the muon anomalous magnetic moment 

a M = (116 591 630 ± 77) x 10~ n . (9) 



We also show that this method is well suited to compute higher-order corrections to new- 
physics contributions to magnetic and electric dipole moments. The only necessary basic as- 
sumption is that the new physics gives a one-loop contribution to the dipole moments, but does 
not significantly affect at tree level composite four-fermion interactions. Under this assump- 
tion, the leading logarithmic two-loop contribution is determined by infrared effects and it can 
be computed from the renormalization-group evolution of the effective theory below the weak 



2 



scale. We find that the QED improvement reduces the one-loop new-physics effect by 6% in the 
case of a M and 11% in the case of the electric dipole moment (EDM) of the electron, assuming 
that the new physics lies around the weak scale. 



2 Two-Loop Calculation of a™ 

The leading two-loop contributions to a^ w come from large QED logarithms. These terms 
correspond to ultraviolet divergences in the effective theory obtained by integrating out the 
heavy modes, and they can be derived in terms of the anomalous dimensions of the dipole and 
current-current operators. 

We start by defining an effective theory valid below the electroweak scale, in which the W 
and Z bosons and the top quark have been integrated out. The effect from heavy particles is 
reflected in higher- dimensional operators. Here we are interested in the operator corresponding 
to a„ 



H » = ~y(£p m » ncjVp ^ Fup ' (10) 

and to other possible dimension-six operators that mix under QED renormalization with H^. 
Since QED is parity conserving and is parity even, it is easy to realize that we need to 
consider only the following parity-even four-fermion operators 



1 1 



Vrf = w"n fluf , Atf = [i^^ii fjulsf ■ (12) 

In eq. ([12]) / indicates a generic fermion different from \x and the factor 1/2 in the definition 
of the operators and compensates the symmetry factor for two identical currents in the 
Feynman rule. The relevant part of the effective Lagrangian is 

C eS = -2V2G ll Y^C i O i , (13) 

i 

where the operators Oi are given in eqs. (PH)-flT^) (with / = {b, r, c, s, d, u, e}) and the Wilson 
coefficients Cj at the electroweak scale \i — M are 



l + J(l-44) 2 



(14) 



CvlM) 



C V JM) = --(l-4s 2 w )(T f - 2Q f s 2 w 



> , C A JM) 



(15) 
(16) 



In eq. (|16|) Tf and Qf are the third isospin component and electric charge of the fermion /, 
respectively. The coefficient Ch is obtained by one-loop integration, while the coefficients in 
eqs. (15)— (16) correspond to tree-level Z exchange. The W boson can only generate operators 
with / = v which are irrelevant for our analysis, because they cannot mix under QED with 
as neutrinos carry no electric charge. 



We are interested in the Wilson coefficient of the operator H„ at the scale /i 



m 



The 



leading-log evolution of the Wilson coefficients from the scale M to a generic scale \i is given 
in terms of the one-loop anomalous dimension matrix 7 and the one-loop beta function. In our 

case, 



e O) , 7 T (e) 
exp / de — — 

e(M) /3(e) 




2b 



V 



cam) 



C^M) , 



(17) 



where e = y/4ira is the QED gauge coupling, and the rotation matrix V is defined such that 
7 = V~ lr y T V is diagonal. The beta function is (3(e) = — b e 3 / (16ir 2 ) , with the coefficient b given 
by 



f 



In eq. (0) the sum is extended over all fermions / with mass less than the scale /1, electric 
charge Qf and multiplicity Nf (Nj = 3 for quarks and Nf = 1 for leptons). Expanding eq. (|T7|) 
in powers of a, we obtain 



ChM = C H ^(M) -^(O^Hj^P-ln-CoAM) , 



(19) 



where ^(Oi^H^) = ^OiH^- Aside from an overall factor G M m^/(2\/2vr 2 ), the first term in the 
r.h.s. of eq. ([19|) gives a^ w (lloop), while the second term gives a^ w (2 loop)^. 

The last ingredient necessary to complete the analysis is the computation of the elements 
^(OijH/j) of the anomalous dimension matrix. This requires a calculation of the divergent 
parts of loop diagrams generating H^, in which a single operator Oi is inserted and a single 



photon is exchanged. Actually, this calculation is completely analogous to the one of the QCD 
anomalous dimension matrix for the AB = 1 effective Lagrangian, relevant for flavour-violating 
bottom-quark transitions. For such a processes a complete list of the divergent contributions 
from the various diagrams can be found in ref . . All we need to do is the proper translation 
from quarks to muons, and from gluons to photons. In this way we obtain 



16 



(20) 
(21) 
(22) 

1 (V, f ,H,) = fQ)N f (23) 
1 (A, f ,H,)=A8Q}N f . (24) 
These expressions are valid both in dimensional regularization with the 't Hooft-Veltman pre- 



40 
~3~ 

,A rr ^ 808 



scription for 75 [I4"| and in dimensional reduction [Id"]. In both schemes there are no finite 
operator renormalizations from non- vanishing matrix elements. 

Replacing in eq. ([19]) the matching conditions of eqs. (|14T)-(|T6|) and the anomalous- dimension 
elements of eqs. (|2l3|)-(p4|), and choosing M = Mz as high-energy scale, we obtain 

5G M rnlaim^) 



EW 



(21oop) LL 



24V2V 



7T 



+ 



36 



43 

T 
2 

27 



215 



In 



m. 



Tf + ^(T f -2Q f s 2 w 



4s 



w > 



In 



m f 



(25) 



where the sum is extended over the fermions with a mass threshold between Mz and m M , 
F = {b, t, c, s, d, u}. Eq. (|^) confirms the results of refs. ||, [L0|, (TTJ, but it disagrees with the 
one presented in ref. ||. It also shows the complete dependence, extending the results of 
refs. H in which the fermionic contribution has been computed in the limit = 1/4. 



The contribution from light quarks is not appropriately described by eq. (|25|), since per- 
turbation theory is not justified. Nevertheless, here we will parametrize the effect by taking 
in eq. ( ^5l) tuq = m u = rrid = m s = 0.3 GeV, and varying tuq between and 1 GeV. In 
this way, our result is consistent with the estimate of the light-quark contribution given in 
ref. [ 10 1 , based on a chiral effective Lagrangian. Taking m c = 1.5 GeV and mj = 4.5 GeV, 
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we find a^ w (2 1oop)ll = —(37 ± 1) x 10~ n . This corresponds to ~ 19% of the one-loop con- 
tribution. We recall that the result in eq. (p5|) includes only the logarithmic contribution. 
The terms not enhanced by large logarithms have been computed in ref. [[| [l l] and amount 
to a^ w (2 loop)jvz, = —(6 ± 2) x 10~ n , where the central value corresponds to a Higgs mass 
Mh = 150 GeV and the error to a variation of Mh in the range between 100 and 1000 GeV. 



3 Three-Loop Calculation of a 



EW 



In the previous section, we have computed a^ w (2 1oop)ll by expanding eq. (|17|) at the first 



order in a. However eq. (|17D also contains the information of higher-order terms, since it 
resums all leading logarithms. We can therefore easily evaluate the magnitude of higher-order 
corrections. For this purpose it is sufficient to expand eq. ( |T7D to the next order in a, because 
{oi/n) \n.{Mz/m^) is much smaller than one. The corresponding result will give us the leading 
logarithmic part of a^ w (31oop). 

Retaining 0(a 2 ) terms, the expansion of eq. (|T7|) yields 



v— a I . a(a) . M 



In — 



47T 



Cj(M) 



(26) 



Because of the presence of the (77) factor, we now need information on the complete structure 
of the anomalous dimension matrix, and not only on the elements given in eqs. fl2"0|)-(|2~3D. We 
start by defining the basis for the required parity-even operators. Besides the operators defined 
in eqs. (|l0|)-([l2l), we need the following four-fermion operators 



Vf=2 f-r»f 



1 



A f=n f^^f flvlsf 



Vff = Hf f'luf , A fr = /Y75/ f'lvlsf with f^f 
Vqq 1 = T1 V 4 q'lvq , Agg> = qYlsq* q'lvim with q ^ q' . 



(27) 

(28) 
(29) 



In eqs. (|27|)-(p9D /, /' (q, q') represent a generic fermion (quark), and all quark operators are 
defined with the colour indices saturated so that each current is an SU(3)c singlet. The 
operators V qq i and A qq > cannot be written in terms of Vff and Ajf with a Fierz rearrangement, 
because of the different colour-index saturation. 
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The matching conditions for the Wilson coefficients at the scale /x = M are 

C Vf (M) = - 2Q f s 2 w ) 2 , C Af {M) = l -Tj , (30) 

Cv ff ,(M) = \{T f -2Q f s 2 w )(T f -2Q f ,s 2 w ) , C Affl (M) = l -T f T f , , (31) 
Cy (M) = Cj (M) = jA qq/ , (32) 

ii ii 4 

where A qq i — 1 if q and q' belong to the same isospin doublet, and A qq > = otherwise. The 
coefficients for the operators of type V, A (V, A) are determined by tree-level Z (W) exchange 
and we neglect Cabibbo-Kobayashi-Maskawa angles. 

The non-vanishing elements of the anomalous- dimension matrix for the operators in eqs. (g7|)- 
(H) are 

16 

l(Vf, Vj) = -Q}(1 + 2N f ) (33) 
1 (V f ,A f ) = 12Q 2 (34) 



52 



1 (A f ,V f ) = -Qj (35) 
16 

l(V fr , V, r ) = -j(Q 2 f N f + Q 2 f ,N f ,) (36) 

j(V ff ,,A ff ,) = j(A ff ,V ff ) = 12 Q f Q f , (37) 

l(V f ,V fr ) = ^Q f Q f (l + 2N f ) (38) 

l{A,,V f f,) = lQ f Q r (39) 
32 

7(y ff ,,Vf) = -QfQ f ,Nf, (40) 
16 

7(%', %») = —Qf'Qf«N f , (41) 

7(VV, V) = 7(4w> %<!') = 12 Wtf (42) 

l(V qql , V q ) = 7 ( V 9 ) = ^QqQq' (43) 

7( V, VW) = 70V, V) = l(Ql + Ql) (44) 

~ ~ 8 

j(Vqq>, Vqf") = l(Aq'i V lf") = ^Qq'Qf" > ( 4 5) 



with/' ^/, g'^and f'^/,/'. 
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Inserting in eq. Q213D the Wilson coefficients at the scale // = M z (see eqs. (P^)-(P^) and 
eqs. (iO) _ dH))' anc ^ the relevant elements of the anomalous-dimension matrix (see eqs. (|20|) - 
(HD and eqs. (H)-©), we obtain 



a* w (31oop) iL 



5 G u m 2 



aim 



1 2 



2 l v '2 " 

Here A and B come respectively from the (77) /2 and the 67 term in eq. 

m u = ma — m Qi are given by 
M z 7826 



(46) 



and, in the 



approximation 
A -- 



2827, 2 

hr 



1/4 and m s 
M z 298 



90 
24 



m. 



In" 

45 m T 



3645 



In" 



M, 
m b 



+ 



7040 , 9 M z 2108 , 9 M 



729 



In 



m n 



+ 



405 



In 2 



m Q 



+ - 



, M z M z 72 M z M z 96, M z M z 
In In 1 m In In m 



m. 



m T 



m, 



rrir 



m,. 



48, M z M z 
— In In 

5 m Q 



128 M z M z 
In In 



1215 



m b 



m. 



(47) 



B 



179 /l 



+ 



45 \3 m b 
2 /, , wij, 4 
77V 3 



+ ln 2 



m T 



+ - In 



M, 



5 
6 

5 V3 



l n 2!^ + I ln 2!^ + 21n 2^ +21n 2 



m c 



3 m c 



2 In 2 

m b 



M2 
m Q 



+ 2 In' 



M, 



m, 



8 / , 9 m c ,9m 
- 2 In 2 — + 2 In 2 — 

5 y mg m 



6/4,9 m T , 9 m T , 9 m T 
+ - - In 2 — + 2 In 2 — + 2 In 2 — 



m. 



m Q 



m, 



8 ln2 mQ 
5 



m,, 



(48) 



Using the same values for the quark masses as in Sect. 2, we find 

a^ w (31oop) LL „„« , M* 



aEW(21oop) Li 



-0.8 ^ In ^ 

7r m M 



which corresponds to a 1% reduction of a® (2 loop) and gives (31oop)x,£ 
Including all the different contributions, we obtain 

a™ = (153 ± 3) x 10- 11 . 



(49) 

0.5 x 10 -11 . 
(50) 



4 New Physics Effects in a^ w and in the Electron EDM 

The effective-Lagrangian method is well suited to discuss effects from new physics. Indeed, in 
the presence of new interactions, characterized by a mass scale A^p of the order of the weak 
scale or larger, only the ultraviolet behaviour of the theory is modified, while the infrared one 



S 



is unaffected. This means that information about new physics can be completely included in 
the matching conditions of the Wilson coefficients, while the QED renormalization described 
by the anomalous dimension matrix through eq. flT7f) remains the same. 

Furthermore, it is reasonable to assume that new physics has a chance to affect sizably 
only matching conditions of operators that, in the SM, are not generated at the tree level, but 
possibly at the quantum level. This is indeed what happens in most of the SM extensions 
generally considered. In this case, the analysis is particularly simple and predictive. 

Let us first consider new physics contributions to a^. Following our assumption, we expect 
that only the matching condition of the operator is affected, but not those of the four- 
fermion operators V and A. This means that the total a M is just given by the sum of the SM 
result discussed in the previous section and of a new contribution a^ p . If a^ p is known at one 
loop, the leading contribution at two loops, enhanced by large QED logarithms, can be simply 
obtained 

< p (21oop) = -l^ln^ af (Hoop) . (51) 
^ 7T m M ^ 

The coefficient in eq. ([5l]) corresponds to the anomalous dimension element 7(iJ At ,iJ M ). For 
instance, for Ajyp = 100 GeV, the inclusion of the two-loop contribution reduces the one- 
loop result by 6%. This result is quite model independent, and it can be applied to specific 
models in which a^ p (lloop) is known, like in the case of supersymmetry [16|, light-gravitino 
interactions |17|], compositeness leptoquarks |]T9 1 , and light non- minimal Higgs bosons [20 |. 

Another case in which the QED logarithms turn out to be large is represented by the EDM 
of the electron, d e . In the SM d e is negligible, but in some of its extensions with new sources 
of CP violation, it can lie just below the present experimental upper bound. For instance, this 
can be the case in some versions of the supersymmetric model |21 . 



The analysis of the QED renormalization can be done along the same lines followed for a M . 
The electron EDM corresponds to an operator in the effective Lagrangian 

- % -d e ea vplb eF v ". (52) 

It is easy to show that in the effective theory below the weak scale, its one-loop anomalous 
dimension is equal to that of the magnetic dipole operator that is given in eq. (|20|) . Therefore 
if d^ 1 describes the new physics contribution obtained by integrating out the heavy modes with 
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mass A^p, the QED improved result is 



4a A NP 



7r m e 

For A NP = 100 GeV, the term inside brackets amounts to a reduction of df^ of 11%. 



(53) 



The result in eq. ([53]) is rather model independent. However, it should be remarked that it 
is valid only if A^p is not much larger than the weak scale. An analysis at very large scales 
should include the mixing of the operator in eq. (|5^) with other CP violating operators, like 
the analogue of eq. (|5*2] ) for the different electroweak gauge bosons, and the analogue of the 
Weinberg operator |22[] for the SU(2) gauge theory. The corresponding anomalous dimension 
matrix can be extracted from ref. |23|]. The result however is more involved than the one 
presented in eq. (§3]), since it depends on the separate unknown coefficients of the different 
operators. 

We wish to thank G. Altarelli and M. Ciuchini for useful discussions. 
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